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We compute the matrix elements of the electromagnetic (EM) operator sF^a^" d between kaon 
and pion states, using lattice QCD with maximally twisted-mass fermions and two flavors of dynam- 
ical quarks (Nf — 2). The operator is renormalized non-perturbatively in the RI'/MOM scheme 
and our simulations cover pion masses as light as 270 MeV and three values of the lattice spac- 
ing from ~ 0.07 up to ~ 0.1 fm. At the physical point our result for the corresponding tensor 
form factor at zero-momentum transfer is /^(O) = 0.417(14 st at)(5 sys t), where the systematic error 
does not include the effect of quenching the strange and charm quarks. Our result differs signifi- 
cantly from the old quenched result /^"(O) = 0.78(6) obtained by the SPQ c dR Collaboration with 
pion masses above 500 MeV. We investigate the source of this difference and conclude that it is 
mainly related to the chiral extrapolation. We also study the tensor charge of the pion and obtain 
the value fx n (0) — 0.195(8 s t a t)(6 sys t) in good agreement with, but more accurate than the result 
/t^O) = 0.216(34) obtained by the QCDSF Collaboration using higher pion masses. 

PACS numbers: 11.15.Ha,12.38.Gc, 13.20.Eb 



INTRODUCTION 



tensor form factor at zero-momentum transfer is 



f**(0) = 0.417 (14 stat ) (5 syst ) = 0.417 (15) 



Accurate measurements of hadron weak decays can 
constrain the parameters of the Standard Model (SM) 
and can place bounds on New Physics (NP) models. In 
particular, the rare decays of the kaon are ideally suited 
to search for new, possibly large CP-violating effects in 
the light-quark sector (see Ref. [![). 

In this work we present a lattice study of the matrix 
elements of the the electromagnetic (EM) operator be- 
tween kaon and pion states, which may be relevant in 
the CP-violating part of the K —> tt£ + £~ semileptonic 
decays. The study has been performed using the gauge 
configurations generated [2} by the European Twisted 
Mass Collaboration (ETMC) with Nf — 2 maximally 
twisted-mass fermions 0, 0| and preliminary results have 
been presented already in Ref. [5J . 

The EM operator involved in the weak s — > d transi- 
tion is given by s F^a^d, where is the EM field 
tensor. Therefore its matrix elements between kaon and 
pion states involve the ones of the weak tensor current, 
which can be written in terms of a single form factor, 



/f% 2 ), as 



(7T \sa- 



flV 



d\K°) = {vWk - «) 



M K + Mir 



(1) 



(2) 



where the systematic error does not include any estimate 
of the effect of quenching the strange and charm quarks. 
Our finding differs significantly from the old quenched re- 
sult fg*(Q) = 0.78(6) obtained in Ref. || by the SPQ cd R 
Collaboration with pion masses above ~ 500 MeV. The 
reason is mainly due to the non-analytic behavior of the 
tensor form factor /^(O) in terms of the quark masses 
introduced by the factor (M^ + M T ) _1 in the parameter- 
ization J!} . Such a behavior was not taken into account 
in Ref. [6j (see later on). 

In the case of the degenerate n — > ir transition, us- 
ing the predictions of the Chiral Perturbation Theory 
(ChPT) carried out in Ref. Q, we obtain for the tensor 
form factor f^(0), known as the tensor charge of the 
pion, the following value 



/r(o) 



0.195 (8 

slat. ) (6 syst ) = 0.195 (10) 



(3) 



which improves the result fj^(0) = 0.216(34) obtained 
by the QCDSF Collaboration [8| with simulations at 
higher pion masses. 



K -> 7T RESULTS 



where q = (px — p-n) is the 4-momentum transfer and the 
factor (Mjf + M^) -1 is conventionally inserted in order 
to make the tensor form factor dimensionless. 

Our simulations cover pion masses as light as 270 MeV 
and three values of the lattice spacing from ~ 0.07 up to 
~ 0.1 fm. At the physical point our result for the K ir 



We have performed the calculations of the relevant 2- 
point and 3-point correlation functions using the ETMC 
gauge configurations with Nf — 2 dynamical twisted- 
mass quarks generated at three values of the lattice 
coupling [3, namely the ensembles A2 — A4 at j3 = 3.8 
(a = 0.098(4) fm), Si - B 7 at /3 = 3.9 (a = 0.085(3) 
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fm), and C x - C 3 at /3 = 4.05 (a = 0.067(2) fm). The 
pion mass ranges between ~ 270 MeV and ~ 600 
MeV and the size L of our lattices guarantees that M V L 
is larger than ~ 3.3. For each value of the pion mass 
and of the lattice spacing we have used several values of 
the (bare) strange quark mass m s to allow for a smooth, 
local interpolation of our results to the physical value of 
m s (see Ref. [9]). The calculation of the 2- and 3-point 
correlation functions has been carried out using all-to- 
all quark propagators evaluated with the one-end-trick 
stochastic procedure and adopting non-periodic bound- 
ary conditions which make arbitrarily small momenta ac- 
cessible. All the necessary formulae can be easily inferred 
from Ref. [Tcj |. where the degenerate case of the vector 
pion form factor is illustrated in details. For each pion 
mass the statistical errors are evaluated with the jack- 
knife procedure. 

The tensor current was renormalized non- 
perturbatively in the RF/MOM scheme in Ref. llj . 



including 0(a 2 g 2 ) corrections evaluated at one-loop in 
lattice perturbation theory {12J. The numerical values 
used in our analyses for the tensor renormalization 
constant are Z T (MS, 2 GeV) = 0.733(9), 0.743(5), 
0.777(6) for /3 = 3.8, 3.9, 4.05, respectively. 

At each pion and kaon masses we determine the tensor 
form factor f^ v {q 2 ) for several values of q 2 < q 2 nax — 
(Mk — M-x) 2 in order to interpolate at q 2 — 1 . Note 
that, because of the vanishing of the Lorentz structure in 
Eq. (p}, it is not possible to determine fr^iq 2 ) at q 2 — 
Imax- We take advantage of the non-periodic boundary 
conditions to reach values of q 2 quite close to q 2 =0. The 
momentum dependence of f^ v (q 2 ) can be nicely fitted 
either by a pole behavior 



Kt, 



(4) 



It V) = /t t (0)/(1 - 
or by a quadratic fit in q 2 

/tV) = /t'(0) • (1 + 4" I 2 + 4" ? 4 ) • (5) 

The good quality of both fits is illustrated in Fig.[TJ where 
the results obtained at two different lattice volumes are 
also compared. It can clearly be seen that: i) finite 
size effects are well below the statistical precision of our 
lattice points; ii) the results for /^(O) and (to a less 
extent) for the slope sf", obtained using the pole dom- 
inance l|4"]). differ only slightly from those corresponding 
to the quadratic fit §5§ in q 2 . 

The values obtained for /^ 7r (0) and s^" 7r depend on 
both pion and kaon masses. The dependence on the latter 
is shown in Fig. [2] for ffi*(0) at AU ^ 450 MeV and it 



1 We remind that the tensor form factor fj$ v (q 2 ) computed in this 
work is C?(a)-improved thanks to the use of maximally twisted 
Wilson quarks 
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Fig. 1: The tensor form factor fx^iq 2 ) obtained at M v ~ 310 
MeV and Mk — 550 MeV versus q 2 in physical units. The 
dots and the squares (shifted for better clarity) correspond to 
the gauge ensembles B\ and Bi , respectively, which differs 
only for the lattice size. The solid and dashed lines are the 
results of the fits based on Eqs. ^ and fBj), respectively. 



appears to be quite smooth. Thus an interpolation at 
the physical strange quark mass can be easily performed 
using quadratic splines. This is obtained by fixing the 
combination (2Mj c — M 2 ) at its physical value, obtaining 
at each pion mass a reference kaon mass, M T ^ ' , given by 



2[M r K ef } 2 - M 2 = 2[M* S ' 2 



{ M phy 



si 2 



(6) 



with M?" 8 = 135.0 MeV and M 



f phys 
K 



494.4 MeV. 
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Fig. 2: Results for fr w (0) versus M 2 K at M n ~ 450 MeV, ob- 
tained assuming the pole dominance ^ for describing the q 2 - 
dependence of our data. The square corresponds to the value 
of fx (0) obtained by local interpolation based on quadratic 



splines (dotted line) at the reference kaon mass M 
MeV from Eq. jgp. 
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The results for /^(O) and sf 1 , interpolated at the 
reference kaon mass M r ^ , are collected in Table [1] and 
shown in Figs. [3] and 0] respectively, for the three lattice 
spacings of our simulations. It can clearly be seen that 
discretization effects are sub-dominant. 
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B 4 (3.9, 24) 


490 


596 


0.681 (11) 


0.671 (11) 


B 5 (3.9,24) 


600 


645 


0.751 (19) 


0.736 (17) 


B 6 (3.9, 32) 


272 


522 


0.537 (16) 


0.529 (16) 


B 7 (3.9,32) 


311 


533 


0.542 (17) 


0.538 (16) 


C*i(4.05,32) 


302 


530 


0.545 (18) 


0.539 (19) 


C* 2 (4.05,32) 


417 


568 


0.648 (17) 


0.649 (16) 


C 3 (4.05, 32) 


486 


594 


0.702 (18) 


0.687 (17) 



Table 1: Results for fr^JO) , obtained at the three lattice spac- 
ings of our simulations using the pole or quadratic £5p 
fits and interpolated at the reference kaon mass for each 
simulated pion mass. 
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Fig. 3: Results for /t*"(0) versus Ml at M K = M r K ef in 
physical units, assuming the pole dominance for describing 
the q 2 -dependence of our data. The dots, squares and triangles 
are our results for the three lattice spacmgs of the ETMC 
simulations, specified in the inset. The solid and dashed lines 
correspond to the fit based on Eq. ^ with B — and D = 0, 
respectively. 

In Ref. [6| the first lattice calculation of the EM op- 
erator matrix element between kaon and pion states was 
carried out in the quenched approximation and for pion 
masses above ~ 500 MeV. There the chiral extrapolation 
was performed adopting a simple linear fit of /p' 7r (0) in 
terms of the squared kaon and pion masses, obtaining at 
the physical point the result f£*(0) = 0.78(6). 

In the degenerate case Mk = M n the chiral expan- 
sion of the tensor current has been studied in Ref. Q- 
The main finding is that the form factor f^ v (0) vanishes 
like Mn for M T -» 0, so that the ratio f^(0)/M n tends 
to a non- vanishing value in the chiral limit. The same 
argument is expected to hold as well in the case of the 
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Fig. 4: The same as in Fig. but for the slope of the 
tensor form factor at q 2 = 0, obtained assuming the pole dom- 
inance for describing the q 2 -dependence of our data. 

K — > 7r transition: the form factor /^"(O) must vanish 
in the SU(3) chiral limit in such a way that the ratio 
f*f*(0)/(M K + M r ) remains finite as (M K , M w ) -> 0. 

Therefore we perform the chiral extrapolation of our 
lattice data using the Ansatz 

/f-(0) = (M r K af + AU)A [l + BMl\og{Ml) + 

CMZ + DM*] , (7) 

where A, B, C and D are low-energy constants, depend- 
ing on the strange quark mass, which are kept as free 
parameters in our fits. The results of the fit assum- 
ing either B = (no chiral log) or D = are shown 
in Fig. [3] by the solid and dashed lines, respectively. It 
can be seen that the effects of the chiral log cannot be 
appreciated with our data. 

At the physical point, after averaging the results ob- 
tained with and without the chiral log in the fitting func- 
tion (JT} and assuming for the momentum dependence ei- 
ther the pole (HJ or the quadratic ([5]) functional forms, 
we get 

/f w (0) = 0.417 (14 stat ) (ETMC) , (8) 

where the error is statistical only. 

Had we neg lected the factor (M^ f + M„) in Eq. © 
the result at the physical point would change only 
marginally: /^(O) = 0.446(19 sta t). On the contrary, 
in the case of the quenched data of Ref. [6j , which were 
determined at pion masses above ~ 500 MeV, the inclu- 
sion of the factor (Mk + M w ) in the chiral extrapola- 
tion changes significantly the result at the physical point 
by many standard deviations, namely from /^^(O) = 
0.78(6) to 

/P(0) - 0.48 (4) (SPQ cd R) . (9) 

These findings indicate that the effect of the quenched 
approximation does not exceed 15%, provided the correct 
mass factor is included in the chiral extrapolation. 
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In the case of the slope Sj" no mass factor should be 
considered in the chiral extrapolation of the lattice data 
shown in Fig. 2J After averaging the results obtained 
assuming for the momentum dependence either the pole 
(UJ or the quadratic ([5]) functional forms, we get at the 
physical point the value 



= 1.10 



t) GeV" 



(ETMC) , (10) 



which is consistent with the quenched result s^" = 
1.11(5) GeV- 2 from Ref. @. 

We now present our estimates of the systematic effects. 

Momentum dependence. We have fitted the (in- 
dependence of our data using either the pole ((4)) or the 
quadratic ([5]) fits (see Fig. [TJ. We have found a system- 
atic error due to these different choices equal to 0.002 for 
/f "(0) and 0.09 GcV~ 2 for sfK 

Finite Size. The comparison between the simulations 
corresponding to the gauge ensembles B\ and B7 (see 
Table Q] and Fig. [IJ, which differs only for the lattice 
size, indicates a very small volume effect of the order of 
0.4% on fg*(0) and of 4% on at M„ ~ 310 MeV. 
We quote therefore a systematic error due to finite size 
effects equal to 0.002 for /f ff (0) and 0.05 GeV" 2 for sf T . 

Discretization. The results obtained using only the 
data at a single lattice spacing a ~ 0.085 fm, where the 
number of simulated pion masses ensures the most reli- 
able chiral extrapolation, are /^(O) = 0.420(10 s t a t) and 
4" = 1.12(6 sta t)- We estimate therefore a systematic er- 
ror due to discretization effects equal to 0.003 for f^^(0) 
and 0.02 GeV~ 2 for s$*. 

Chiral extrapolation. As shown in Figs. [3]|4l the differ- 
ence in the extrapolations to the physical point including 
or excluding a chiral log in the fitting function is quite 
small, being equal to 0.002 for fg K (0) and 0.02 GeV~ 2 
for 7r . We take these numbers as our estimates of the 
uncertainty due to the chiral extrapolation. 

Adding all the systematic errors in quadrature, our 
final results are 



/**(0) = 0.417 (14 stat ) (5 syst ) 
= 0.417 (15) , 



Kir 

Of 



1.10 (8 

slat 

) (ll syst ) GeV" 2 
1.10 (14) GeV~ 2 , 



(11) 



(12) 



where the systematic errors do not include any estimate 
of the effect of quenching the strange and charm quarks. 



off-diagonal entries in the down-type squark mass ma- 
trix, already calculated in Ref. [6J using the quenched 
estimate /f^O) = 0.78(6). 

The master formula, expressing the CP violating part 
of the contribution of the EM operator to the rare decay 
Kl — > 7r°e + e~ , can be written as (see Ref. @) 



B(K L -> 7r u e + e-) = 5.3 ■ 10" 4 B T (Im 8+Y 



(13) 



where the numerical coefficient (5.3 • 10~ 4 ) is the ap- 
propriate one for the EM operator renormalized in the 
MS scheme at a renormalization scale of 2 GeV and 
for a gluino and average squark masses of 500 GeV. In 
Eq. ()13[) the quantity Bt is the main hadronic parame- 
ter related to the tensor form factor at zero-momentum 
transfer /f "(0) by 



2M 



K 



/t*(0) 



■ T 



(14) 



where f+*(0) is the vector form factor at zero- 
momentum transfer, appearing in the semileptonic 
decay K + — > 7r°e + ^ e , while T is a correction which takes 
into account the different g 2 -dependence of the vector 
and the tensor form factors (see Ref. [6j for its defini- 
tion). The quantity T can be related to the slopes s+*' 
and sf T of the vector and tensor form factors at q 2 = 



by the following relation 

T = 1 — 2M 2 (s^ 



R 



(15) 



where R is a simple kinematical factor depending on the 
kaon and pion masses. At the physical point one has 
R = 1.855. 



In Ref. [14[ we investigated the vector and scalar form 
factors entering the Km semileptonic decay, obtaining 
the results f^(0) = 0.9560(57 sta t)(62 syst ) and s 1 ^ = 
1.30(13 s tat)(12 sy st) GeV -2 . Using the latter values and 
the results (jlH12|l . the correction factor T turns out to 
be quite close to unity, namely T — 0.986(10 st at)(llsyst), 
and for the quantity Bt we get 



B T = 0.676 (24 stat ) (12 syst ) = 0.676 (27) 



(16) 



Therefore, using in Eq. (|13[) the experimental upper 



bound B(K L ->■ 7r°e+e-) < 2.8 • 10 



-10 

and our result (|16p . we obtain the bound 

|Im S + \ < 1.1- 1Q- 3 (90% C.L.) 



C.L.) [131 
(17) 



Bound on the supersymmetric coupling 5+ 

We want now to use our new determination (fTT|) 
for the tensor form factor f^ jO) as well as the re- 
cent experimental bound [13j on the branching ratio 
B(Kl —> 7T°e + e~) to update the bound on the super- 
symmetric coupling related to the splitting of the 



PION TENSOR CHARGE 

Following Ref. [7| the chiral expansion of the pion ten- 
sor charge /J*'(0) has the form 



fF(0) - M. A' r , (4?r/)2 
C'Ml+D'M^} , 



M 2 

7^ log(M 2 ) + 



(18) 
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where / ~ 122 MeV [9j is the pion decay constant in 
the SU(2) chiral limit and the presence of the factor M v 
is expected to have an important, bending effect on the 
value extrapolated to the physical point. 

Our results for f^(0), obtained at three values of the 
lattice spacing and for 270 MeV <M T < 600 MeV, are 
shown in Fig. [5] and compared with the Nf = 2 results of 
Ref. [I], having > 440 MeV, and with the quenched 
calculations of Ref. 6] , ranging from M„ ~ 530 MeV up 
to Mn ~ 800 MeV. It can be seen that our results have a 
much better statistical precision and cover much lighter 
pion masses, where the bending effect due to the overall 
factor Mtt is clearly visible. 



Ensemble 


M n 
(MeV) 


/T(0) 
(pole) 


/T(0) 
(quadratic) 


A 2 (3.8,24) 


422 


0.568 (15) 


0.555 (13) 


As (3.8, 24) 


491 


0.643 (17) 


0.623 (16) 


A 4 (3.8,24) 


598 


0.780 (18) 


0.767 (17) 


Bi(3.9,24) 


319 


0.448 (16) 


0.432 (14) 


B 2 (3.9,24) 


393 


0.500 (18) 


0.487 (18) 


B 3 (3.9,24) 


453 


0.586 (22) 


0.580 (20) 


B 4 (3.9,24) 


490 


0.661 (24) 


0.641 (22) 


B 5 (3.9,24) 


600 


0.736 (19) 


0.720 (18) 


B 6 (3.9,32) 


272 


0.384 (23) 


0.371 (20) 


B 7 (3.9,32) 


311 


0.431 (18) 


0.426 (15) 


Ci (4.05, 32) 


302 


0.420 (26) 


0.404 (23) 


C* 2 (4.05,32) 


417 


0.569 (16) 


0.562 (15) 


C 3 (4.05, 32) 


486 


0.654 (19) 


0.640 (17) 



I, 



Table 2: 

ings of our simulations IQj, using the pole or quadratic fits for 
describing the q 2 -dependence of our data. 
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Fig. 5: Results for the pion tensor charge /^^(O) versus M 2 
in physical units from our simulations at the three different 
values of the lattice spacing (open markers), obtained assum- 
ing the pole dominance ^ for describing the q 2 -dependence 
of our data, and from Refs. 0/ (full triangles) and Jjfy (full 
squares). The dashed line is the result of the fit H18V with 
D' 7^ applied to all ETMC points, while the solid line cor- 
responds to the use of the fitting function 118}) with D' = 
applied to the ETMC data with M n < 420 MeV. The dotted 
line represents the fit described in the text and applied to the 
quenched data of Ref. Q/. 



final result is 

fi"(0) = 0.195 (8 stat ) (6 syst ) 

= 0.195 (10) (ETMC) (19) 

which improves the QCDSF result Q 

/^(O) =0.216 (34) (QCDSF). (20) 



Applying Eq. (Tig)) with D' = to our lattice points 
with M„ < 420 MeV (see the solid line in Fig.©, we get 
at the physical point f$*(0) = 0.195(8 sta t)- 

As for the estimates of the systematic effects, we follow 
the same procedures described in the previous Section for 
the K — > 7T transition, obtaining an effect equal to 0.004 
from fitting the momentum dependence with either a pole 
or a quadratic forms, 0.003 from the finite size 2 , 0.001 
from discretization and 0.003 from the chiral extrapola- 
tion. Adding all the systematic errors in quadrature, our 



2 The comparison between the simulations corresponding o the 
gauge ensembles B\ and B7, which differs only for the lattice 
size, indicates a volume effect of the order of 2 — 4% on /^(O) 
(see Table [2}, which is close to the statistical error. Such an 
effect is however mostly due to the volume effect on the pion 
mass. Indeed, the finite size effect on the ratio f^, 7 ' (0) / 'M n re- 
duces to ~ 1.5%, which we therefore take as our estimate of the 
corresponding systematic error on our final result. 



Finally, we also apply a simple fit of the form /^(O) = 
M v A' [l + C'Af 2 ] to the three quenched data of Ref. @ 
(see the dotted line in Fig. [5]), obtaining at the physical 
point the result 

/? ff (0) = 0.221 (21 stat ) (SPQ cd R) , (21) 

which clearly shows that quenching effects on the pion 
tensor charge are sub-dominant with respect to the 
present precision. 

Before closing we want to comment briefly on the 
SU (3) breaking effect on the tensor form factor at 
zero-momentum transfer. Our final values obtained for 
/f ""(0) and /^(0) [i.e., Eqs. (n]) and (19])] differ by a 
factor of w 2. This huge difference is however mostly 
due to the factor (Mk + M n ), conventionally inserted in 
Eq. |l|. By introducing the quantities 

F^(0) = f«*(0)/(M K + M n ) (22) 
F$*(0) = fF{Q)/2M„ (23) 



6 



one gets from Eqs. (jlll) and (IT9|) the values 

F^(0) = 0.663 (22 stat ) (8 syst ) GeV" 1 

= 0.663 (24) GeV" 1 (24) 

F^(0) = 0.722 (30 stat ) (22 syst ) GeV" 1 

= 0.722 (37) GeV" 1 (25) 

which yield a ratio F^ 7T (0)/Ff n (0) equal to 0.92(6). This 
indicates that the SU(3) breaking effect on the matrix el- 
ement of the EM operator sF^a^d at zero-momentum 
transfer is equal to —8(6)%. This value is of the same 
sign and order of magnitude of the SU (3) breaking ef- 
fect on the corresponding matrix clement of the vec- 
tor current, namely on the vector form factor at zero- 
momentum transfer f^ v (0). Indeed, taking into account 
that f+*(0) = 1 because of charge conservation and 
f£ n (0) = 0.9560(84) from Ref. 44|, the SU(3) break- 
ing effect on the matrix element of the vector current at 
zero-momentum transfer is equal to —4(1)%. 
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